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We derive a new QCD sum rule for -D(O^) which has only the Dtt continuum with a resonance 
in the hadron side, using the assumption similar to that has been successfully used in our previous 
work to the mass of Ds(0+)(2317). For the value of the pole mass Mc = 1.38 GeV as used in the 
Ds{Q'^) case we find that the mass of D{Q'^) derived from this sum rule is significantly lower than 
that derived from the sum rule with the pole approximation. Our result is in agreement with the 
experimental dada from Belle. 
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I. INTRODUCTION 

BaBar Collaboration discovered in 2003 a positive- 
parity scalar charm strange meson Dsj(2317) with a very 
narrow width [H, which was confirmed by CLEO later 
In the same experiment CLEO ,2] observed the 1+ 
partner state at 2460 MeV. Since these two states lie be- 
low DK and D*K threshold respectively, the potentially 
dominant s-wave decay modes Z3sj(2317) DgK etc are 
kinematically forbidden. Thus the radiative decays and 
isospin-violating strong decays become dominant decay 
modes. Therefore they are very narrow. 

The discovery of these two states has inspired great 
interest in their nature in literature. The key point is 
to understand their low masses. The (2317) mass is 
significantly lower than the expected values in the range 
of 2.4 — 2.6 GeV in quark models The model using 
the heavy quark mass expansion of the relativistic Bethe- 
Salpeter equation predicted a lower value 2.369 GeV 
which is still 50 MeV higher than the experimental data. 

Later Belle [H observed the wide -0(0"'") resonance with 
mass Mo = 2.308±0.0017±0.0015±0.0028Gev and width 
around F — 276 MeV. Another puzzle arises from this 
result: why are these two resonances nearly degenerate 
in mass while Ds{0~) is 100 MeV higher than D{0-)7 

The earlier results for the mass of Ds{0'^) from QCD 
sum rules are either significantly larger than the exper- 
imentally observed mass of £'sj(2317) [Q] or consistent 
with it within theoretical uncertainty but with signifi- 
cantly larger central value The results from lattice 
QCD calculations are similar, see @ and [1, H^l- 

Recently, there have been two investigations on this 
problem using sum rules in full QCD including the 0{as) 
corrections. In Ref. [Tlj the value of the charm quark 
pole mass = 1.46GeV is used and 0~^cs is found to be 
100 — 200 MeV higher than the experimental data. On 
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the other hand, in Ref. Mc ~ 1.33GeV is used and 
the central value of the results for the 0"*" cs mass is in 
good agreement with the data. As commented by the 
author, the uncertainty in the value of Mc is large. 

The difficulty with the cs interpretation leads many 
authors to speculate that Dsj(2317) is a cqsq four quark 
state [ll,[ij| or a strong Dtt atom [l^. However, quark 
model calculations show that the mass of the four quark 
state is much larger than the 0+ cs state [H, [13 • Fur- 
thermore, the four quark system has two 0"*" states. Only 
one has been found below 2.8 GeV in the experimental 
search, consistent with the cs interpretation. 

It was suggested in Ref. [l8| that the low mass of 
1357(2317) could arise from the mixing between the 0"'" 
cs state and the DK continuum. It was pointed out in 
Q that in the formalism of QCD sum rules the physics of 
mixing with DK continuum resides in the contribution 
of DK continuum in the sum rule and including this 
contribution should render the mass of Ds(0'^) lower. 

Usually the contribution of the two-particle contin- 
uum is neglected in the QCD sum rules. In a recent 
work [1^ we argued that because of the large s wave 
coupling of Ds{0~^)DK [l^, and the adjacency of 
the Ds{0~^) mass to the threshold this contribution may 
not be neglected. We calculated this contribution under 
the assumption that the form factor for the coupling of 
the scalar current to the two-particle continuum in the 
low energy region is dominated by the bubble diagrams 
formed by the coupling of the 0"'" to the two particles. We 
found that including this term in the sum rule indeed ren- 
ders the mass and decay constant of Ds{0'^) significantly 
lower. Using the pole mass value Mc = 1.38 GeV the 
mass of Ds{0~^) is found to be in good agreement with 
experimental data. 

In the present work we would like to show that the 
same assumption for the form factor for the coupling 
of the scalar current to the two-particle continuum also 
leads to the mass value of -0(0"*") in agreement to the ex- 
perimental data, thus, explains the nearly degeneracy of 
these two states. 

In Section|TT]we derive the Dn continuum contribution 
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and write down the full QCD sum rule for D{0~^). The 
numerical results of the mass of -0(0"*") are collected in 
Section Uni 



II. THE QCD SUM RULE FOR THE SCALAR 
CHARM MESON 



We consider the scalar correlation function 
U{q^)=i I d^xexp[iga;](0|T{J(a;), jt(o)}|0) (1) 



where J(x) = c(x)d(x) is the interpolating current for the 
charged scalar charm meson, n(g^) satisfies the following 
dispersion relation 



ds' 



Imn(s) 
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The leading terms of the imaginary part Imn(s) at the 
quark gluon level and its as correction have been calcu- 
lated in [13, • After making the Borel transformation 
to suppress the contribution of higher excited states and 
invoking the quark-hadron duality, one arrives at the sum 
rule 
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where nic and mj^ are the charm and down quark current 
mass. The charm quark pole mass is defined as (23 | 



as 

M2 TT 



(dd) is the down quark condensate, {agC^) is the gluon 
condensate, {ggda ■ Gd) is the quark gluon mixed con- 
densate. Mb is the Borel mass. The radiative correction 
function reads 

G{x) = - + 2Li2(.T) + log(x) log(l -x)-l \og{- - 1) 
4 2 X 



■ log{x) 



- log(l ~ x) + x log{- - 1) - — 

X 1 



The above equation is of the same form as that for 
Ds{0^)- The difference between the two cases is the fol- 
lowing. Dg{0'^) is a very narrow resonance below the 
DK threshold. Neglecting the small iso-spin violating 
interaction it can be represented by a pole term in the 



sum rule. On the other hand, £'(0"'") is a very broad res- 
onance in the Dtt spectra and there is no pole term in 
the hadron spectra. In the traditional treatment of QCD 
sum rules the very broad resonance is still approximated 
by a pole term solely for simplicity and without justifi- 
cation. Moreover, the remaining two-particle continuum 
is completely neglected. 

According to this procedure the spectral density at the 
hadronic level is taken to be the pole term plus the con- 
tinuum starting from some threshold which is identified 
with the QCD continuum. 



ImnH(t) 
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where /o is the vector current decay constant of 0+ cd 
particle analogous to = 132 MeV, Mq is the mass of 
this particle, and sq is the continuum threshold above 
which the hadronic spectral density is modelled by that 
at the quark gluon level. The recent work [l^ also uses 
the above ansatz. The results from the above sum rule 
with Mc = 1.38 GeV are shown in Fig. 1. This may 
not be a good approximation. The importance of Dtt 
continuum in the D(O^) channel was first emphasized 
by Blok, Shifman and Uraltscv in Ref. [22] from the 
consideration of duality. 




FIG. 1: Mo vs Mg from the sum rule with pole approximation 
for D{0'^). Curves from top to bottom correspond to so = 
8.0, 7.5, 7.0 Gev^ respectively. 



In our previous work [131 it was pointed out that the 
strong s wave coupling of the Ds{0'^) with the two parti- 
cle DK state and the adjacency of the Ds{0~^) mass to the 
DK continuum threshold result in large coupling channel 
effect which corresponds to the configuration mixing in 
the formalism of the quark model. Therefore, the two- 
particle continuum can not be neglected. In the present 
case of -0(0"'") the two-particle term is the whole hadron 
spectra needed in the sum rule. We calculate this term 
with the same assumption used in [l^. Besides giving a 
more accurate sum rule for the mass of D(0"'"), this also 
constitutes a check for the approximation used in for 

The contribution of the Dtt continuum to the hadronic 
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spectral density reads 
ImnH(t) 
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(Mq- 
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x\F{t)\''e{Vi ~ Md - m^)e{so - t) 
+QCD continuum x 9(t ~ sq) , (5) 



where F{t) is the form- factor defined by 
F{t) = (0|c(0)d(0)|i:»7r) 



(6) 



Similar to the approximation used in [1^ , from the large 
s wave coupling of D{0'^)D'!t we assume that in low en- 
ergy region F{t) is dominated by the product of a factor 
proportional to the propagator of the D{0~^) and a factor 
from the final state interaction represented by the chain 
of the bubble diagrams generated from the D{0^)Dtt in- 
teraction as shown in Fig. 2. 



its contribution can be neglected. Summing the chain of 
the bubble diagrams, we have 



Fit) 



A 



where 
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From Cutkosky's cutting rule, 

^2 



{q - k)^ - Ml 
-|-subtraction (7) 
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where 



kit) 



^it - iMp + m^)2) it - jMp - m^)2) 
2^/t 



(8) 



and we have taken into account two intermediate states 
of different charge. From the dispersion relation. 



FIG. 2: Heavy, light and dotted line represents -D(O^), D 
and TT respectively. Black circle represents the Born s wave 
amplitude of Dn scattering. Blank circle represents the scalar 
current. 



Let g be the coupling constant in the effective 
£'(0~'")-D7r interaction lagrangian. The Born s wave am- 
plitude in the figure contains both a t channel pole term 

2 

and a crossing s channel pole term 
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Using the relation between s, t and the scattering angle 9 
one finds that the latter is an analytic function of t with 
only a short cut. Therefore, this term can be approxi- 
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mated by a pole form j^zt^, where 
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Different from the case of L's(0+), the crossing term is 
actually small for experimental values of the masses in 
the present case of 15(0"'"). But we shall still keep it in 
the formulae. In the chiral lagrangian the interaction 
between D and tt through exchanging two pions is sup- 
pressed by (iJ/2\/27r/7r)"' for low energy E of tt. Hence 



1 [ ,3g2 fcftO it -Ml) 
n J 1671 {f - MDit' - t) 



where Mo is the renormalized mass of the -D(0+) me- 
son. The mass renormalization condition leads to 
ReT,iMl) = 0. Therefore we have 

Fit) = X{it-Ml-^il+^-^) 



,,, k{t') (t~Afg) N 



xP / 

(MD+m,)2 



(9) 



where P denotes the principal part of integration. 

The unknown constant A in Fit) can be determined by 
applying the soft-pion theorem to the extrapolated value 
of the matrix element (0|c(0)d(0)|i:>7r) at t = M^ with 
the result 



FiMl) 



IdMI 



hirric + md) 
Hence we have 



III. THE MASS OF D(0+) 

In the numerical analysis, we use = 8 MeV, = 
1.18 GeV, A^c = 1-38 GeV, (Jd) = -(0.243)^ GeV^, idg^ 



4 



Gd) = O.SGeV^ * (dd), (a.G^) = 0.06 GeV^ Aqcd = 
0.325 GeV, Md = 1.869 GeV, = 0.140 GeV, fo = 0.2 
GeV = 0.132 GeV [26]. 

The g value has not been determined very well theo- 
retically. It wasj^ound to be in the interval g — 7.5 — 5.1 
GeV in Ref. [20|. Inclusion of the contribution of Dtt 
continuum in the sum rule analysis of the scalar current 
channel may lower the g value. Since the uncertainty is 
large, we determine it from the experimental width F of 
-D(0+) through the equation 

2 _ ^, 3V(Afg - imp + m^ f j) (Mg - (mn - 
^ 32^Af3 J 

and use the experimental central value F = 276 MeV [1| 
as the input. For the experimental central value of the 
D(0+) mass M„ = 2.308 GeV, g = 8.09 GeV. Notice that, 
since we can not do wave function renormalization for 
£'(0"'') at the resonance energy Mq, the coupling constant 
g defined in this work is not exactly the correspondence 
of that defined in Ref. [3 for D,{0+). 
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FIG. 3: The variation of the ratio of the LHS to the RHS of 
the sum rule with Af| when so = S.OGeV^, Mq = 2.297 GeV. 
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FIG. 4: The variation of the ratio of the LHS to the RHS of 
the sum rule with M| when so = 7.5GeV^, Mo = 2.285 GeV. 



The good convergence of the OPE series and domi- 
nance of the Dtt term over the QCD continuum beyond 
So constrain the Borel mass in a region depending on 
Mc and So- For = 1.38 GeV, so = 8.0,7.5,7.0 
GeV2, Ml e [1.15,2.5], [1.15,2.25], [1.15,2.2] GeV^ re- 
spectively. Since the hadron side of the sum rule de- 
pends on the unknown value Mq in a complicated way, 
we have to input a " trial" value of Mq and require that 
the ratio of the hadron side (the left hand side or LHS) 
to the quark-gluon side (the right hand side or RHS) of 
the sum rule is equal to 1 in the middle of the working 
range of M^. The results are shown in Fig. 3, 4, 5 for 
So = 8.0, 7.5, 7.0 respectively. 

Comparing these results to those obtained in the pole 
approximation shown in Fig. 2 one finds that the re- 
sults for the mass of -D(0+) are lower by 80 — 40 MeV. 
For Mc = 1.38 GeV, so = 8.0 - 7.0 GeV^ we find 
Mq = 2.30 — 2.27 GeV in agreement with the experimen- 
tal data. The curve is very sensitive to the trial value 
of Mq, hence the uncertainty of our results from differ- 
ent values of Mb for fixed values of other parameters are 
small. Therefore, with the same approximation for the 
two-particle continuum and same values of Mc and so we 
obtain simultaneously the correct values of the masses 
of the narrow resonance Ds{0~^) below the threshold and 
the broad resonance 15(0"'") above the threshold. 

Under our approximation for the two-particle contin- 
uum and acceptable value of g, there is no solution for 
the sum rule for -0(0"*") containing a below-threshold pole 
term and the two-particle continuum in the hadron side 
as that used in the Ds{Q'^) case. 
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FIG. 5: The variation of the ratio of the LHS to the RHS of 
the sum rule with A4"| when so = 7.0GeV^, Mo = 2.270 GeV. 
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